Constraining crystalline color superconducting quark matter with gravitational-wave 

data 
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We estimate the maximum equatorial ellipticity sustainable by compact stars composed of crys- 
talline color-superconducting quark matter. For the theoretically allowed range of the gap parameter 
A, the maximum ellipticity could be as large as 10 -2 , which is about 4 orders of magnitude larger 
than the tightest upper limit obtained by the recent science runs of the LIGO and GEO600 grav- 
itational wave detectors based on the data from 78 radio pulsars. We point out that the current 
gravitational-wave strain upper limit already has some implications for the gap parameter. In par- 
ticular, the upper limit for the Crab pulsar implies that A is less than O(20) MeV for a range of 
quark chemical potential accessible in compact stars, assuming that the pulsar has a mass 1.4Mq, 
radius 10 km, breaking strain 10 -3 , and that it has the maximum quadrupole deformation it can 
sustain without fracturing. 

PACS numbers: 04.30.Db, 25.75.Nq, 26.60.+C 



Introduction. When nuclear matter is squeezed to a 
sufficiently high density, there is a transition from nu- 
clear matter to quark matter. Since the density required 
for the transition to happen is believed to be not much 
higher than nuclear-matter density, the dense cores of 
compact stars are the most likely places where quark 
matter may occur astrophysically. Except for hot new- 
born compact stars, it is now generally believed that the 
deconfined quark matter (if exists) in the interior of com- 
pact stars is in a color-superconducting phase [H, 0, IE 0] j 
in which the quarks form Cooper pairs due to the BCS 
mechanism (see [5j and references therein for reviews). 
At sufficiently high densities, the favored pairing phase 
is the color-flavor-locked (CFL) phase 'S\, in which pair- 
ing between quarks of different colors and flavors is al- 
lowed. In the intermediate density regime relevant to the 
cores of compact stars, it is found that crystalline color- 
superconducting quark matter is a more favored phase 
[E B H S El- However, it should be noted that so 
far the studies on the crystalline phase are based only on 
phenomenological models of QCD, mainly the Nambu- 
Jona-Lasinio model [ll| . The true ground state of quark 
matter in the cores of compacts stars is still a matter 
of debate. Furthermore, there is as yet no study on the 
construction and stability of a hydrostatic equilibrium 
stellar model composed of crystalline quark core in gen- 
eral relativity. 

While it is still an open question whether quark stars 
exist in nature, it is interesting to ask how their ex- 
istence would affect observable phenomena related to 
compact stars. Furthermore, could we make use of ob- 
servation data to constrain various theoretical models? 
In this paper, we show that the observational upper 
limits on gravitational-wave emission from known pul- 
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sars can be used to constrain compact star models com- 
posed of crystalline color superconducting quark matter. 
The prospect of detecting the gravitational-wave signals 
emitted by quark stars has been considered before (e.g., 
[2 EE EH HE SB])- However, our work represents the 
first attempt to make use of real observational data to set 
constraint on theoretical parameters of the quark matter. 

Rotating compact stars are among the most promis- 
ing gravitational-wave sources for Earth-based interfero- 
metric detectors such as LIGO, VIRGO, GEO600, and 
TAMA300. These objects can emit gravitational waves if 
they are asymmetric about the rotation axis. A number 
of mechanisms to give rise to the asymmetry have been 
proposed, including (1) nonaxisymmetric distortions of 
the solid crust or core of the star [HI EE EE HJ; (2) 
rotational induced instabilities such as the bar-mode and 
r-mode instabilities (see [2l|, [22[ for reviews) ; (3) dis- 
tortion produced by strong magnetic fields 123.1241 HH ; 
and (4) free precession of the star [IE l27l. |28||. Here 
we study the nonaxisymmetric distortions of solid com- 
pact stars composed of crystalline color-superconducting 
quark matter and show that the recent observational re- 
sults by LIGO and GEO600 are already sensitive enough 
to put constraints on the gap parameter A associated 
with the color-superconducting phase. 

Over the past few years, the LIGO and GEO600 de- 
tectors have conducted four science runs starting in 2002. 
In particular, data from the second science run (denoted 
as S2) have been used to set direct upper limits on 
the gravitational-wave amplitudes from 28 known pulsars 
[29l | . The lowest gravitational- wave strain upper limit in 
the S2 run is at the level ~ 10~ 24 and the equatorial 
cllipticities are ~ 10~ 5 . Recently, the LIGO Science Col- 
laboration has presented the latest results for 78 pulsars 
based on data from the third and fourth science runs (de- 
noted as S3 and S4) [3fJ . The improved sensitivity of the 
detectors in the S3/S4 runs leads to the tightest strain 
upper limit ~ 10 -25 and equatorial ellipticity ~ 10 -6 . 
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For conventional neutron star models, the maximum 
ellipticity that can be sup por ted by the neutron star's 
solid crust is about 10~ 7 [ill Il9j |. The current detec- 
tors are still not sensitive enough to probe the upper 
range of ellipticity permitted by neutron stars. How- 
ever, the current limits are well into the range permitted 
by some exotic compact star models. Based on the hy- 
pothesis proposed by Xu [3l| that quarks might form 
clusters via a Van der Waals-type interaction, Owen [l9| 
found that solid compact stars composed of such quark- 
clusters might sustain ellipticity up to ~ 10 -4 . It should 
be noted that the solid quark-star models considered in 
this paper, namely, the crystalline color-superconducting 
quark matter, is more robust than Xu's hypothesis from 
a theoretical point of view. We find that the maximum 
ellipticity sustainable by these models could be as large 
as 5 x 1CP 2 . This makes the current strain upper lim- 
its obtained by the LIGO and GEO600 detectors become 
more astrophysically relevant and interesting. 

Elastic deformations and gravitational-wave emission. 
A triaxial pulsar, rotating about a principal axis, would 
radiate gravitational waves at twice the rotation fre- 
quency. In the quadrupole approximation, the charac- 
teristic strain amplitude is [32[ 



, 16tt 2 GW 2 
ho = 1 e, 
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where / is the pulsar's spin frequency, r is the distance to 
the pulsar, e = (I X x—Iyy)/Izz is the equatorial ellipticity, 
and Iij is the moment of inertia tensor. The z axis is the 
rotation axis. Using the definition of the mass multipole 
moment, Qi m = J pr l Y* m d i x 1 it can be shown that the 



ellipticity is related to the m 
Q22 by 
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In this paper, the pulsar's quadrupole moment is as- 
sumed to be due to its elastic deformation in the presence 
of a solid core. Assuming that the material is everywhere 
strained to the maximum, the maximum quadrupole mo- 
ment the star can sustain is given by Ushomirsky et al. 
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UU + U 2 -r^) dr, 
dr / 



(3) 

where cr max is the breaking strain, v is the shear modu- 
lus, g — GM (r)/r 2 is the local gravitational acceleration, 
and U = 2 + dhig/d\nr. For a uniform-density incom- 
pressible model, which is a good approximation to quark 
stars considered in this paper, <3 max is simplified to 



Qn 



GM ' 



(4) 



It should be noted that in deriving Eq. ([3]), the self- 
gravity of the deformation is neglected. However, drop- 
ping that approximation would only change <3 max by a 



factor of 2 [20|]. This is not large enough to change the 
implications we shall draw from the estimate of the ellip- 
ticity described below. 

It is seen from Eq. ^ that more rigid material (i.e., 
material with a larger shear modulus v) can lead to a 
larger quadrupole deformation, which in turn implies a 
larger ellipticity [see Eq. ([2])]. This is the main reason 
why solid quark stars can sustain much larger ellipticity 
than conventional neutron stars since the shear modulus 
of solid quark matter can be a few orders of magnitude 
larger than that of the neutron star's crust [l|| (also see 
below) . In Eq. © we define the shear modulus to be half 
the ratio of the stress to the strain, which has a factor 
of 2 different from that defined in the original paper by 
Ushomirsky et al. [HI]. This definition is consistent with 
that used by Mannarelli et al. 33] in calculating the 
shear modulus of crystalline color-superconducting quark 
matter, which is given by 
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where A is the gap parameter and [i is the quark chem- 
ical potential. We remark that this result is obtained 
by performing a Ginzburg-Landau expansion to order 
A 2 . Since the control parameter for the expansion is 
[A/(m 2 /8/i)] ~ 1/2 (with m s being the strange quark 
mass) for the favored crystalline phase [33[, Eq. §5§ can 
only be used to fix the order of magnitude of v. Higher- 
order correction terms are suppressed only by about 1/4. 

To describe quark matter within compact stars, the 
relevant range of fx is 350 MeV < fi < 500 MeV [H, B 
In the lower region of this window, Ippolito et al. [34| 
notice that the crystalline phase is more likely to be a 
two- flavor ((ud)), rather than the three- flavor ((ud) and 
(us)) condensate as considered in [33J. In the CFL phase, 
where the density is high enough so that the effects of 
the strange quark mass can be neglected, the gap pa- 
rameter A is expected to lie between 10 and 100 MeV 
UU- However, for the quark matter to be in the 
crystalline phase rather than the CFL phase, Mannarelli 
et al. [HJ estimate that the reasonable range for A is 
5 MeV < A < 25 MeV. These windows for fj, and 
A imply that the shear modulus of crystalline color- 
superconducting quark matter within compact stars is in 
the range 0.47 MeV/fm 3 < v < 24 MeV/fm 3 . For com- 
parison, the shear modulus of the neutron star's crust is 
of the order 1 keV/fm 3 [H,^. 

The breaking strain cr max is highly uncertain. This is 
the value beyond which solid materials no longer behave 
elastically. Instead, they would either crack or undergo 
plastic deformation. For perfect crystals without defects, 
c ma x could be as high as 10 _1 [35|]. However, it is also 
known that c max of real crystals can be several orders of 
magnitude lower. The discrepancy is explained by the 
presence of defects (e.g., dislocations) in real crystals. 
Based on the extrapolation of terrestrial materials, val- 
ues as high as 10 -2 has been suggested for neutron star 
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FIG. 1: (color online). The areas above the solid lines are 
excluded by the direct upper limit for the Crab pulsar ob- 
tained from the S3/S4 runs for the cases <r max = 10~ 3 and 
1CP 2 . The dotted (dashed) line is the constraint set by the 
spin-down limit for er max = 10 -3 (10 -2 ). The rectangular box 
is the theoretically allowed region of fi and A. 



crusts [3 a . |37I | , which may also be favored by the enor- 
mous energy (~ 10 46 erg) liberated in the December 27, 
2004 giant flare of SGR 1806-20 according to the magne- 
tar model [H, [39[ . In this work, we consider the effects 
of Umax in the range 10 -3 to 10~ 2 , values that have been 
used in the study of compact stars by other authors (e.g., 

[H, MM, Hi). 

Using the shear modulus given in Eq. ([5]), the max- 
imum quadrupole moment for a solid quark star in the 
crystalline color-superconducting phase can be estimated 
by Eq. However, the moment of inertia I zz is needed 
to calculate the corresponding maximum equatorial cllip- 
ticity. Using the empirical formula for strange stars given 
by Bejger and Haensel (see Eq. (10) of 40]), Eq. Q gives 
the maximum equatorial ellipticity as 



e max = 2.6 x 10" 
R 
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We use the values M = 1.4M©, R = 10 km, and 
Cmax = 10~ 2 in order to compare with in which 
Owen obtained e max ~ 2 x 10~ 4 for solid strange stars 
(with the quarks clustered in groups of about 18). For 
the estimated range of the shear modulus given above, we 
find that e max could be as large as ~ 5 x 10 -2 for solid 
quark stars in a crystalline color-superconducting phase. 
This relatively large value of e max is about 4 orders of 
magnitude larger than the tightest upper limit obtained 
by the combined S3/S4 result for the pulsar PSR J2124- 
3358 [U. 

Constraints set by the Crab pulsar. Eq. ([T|) suggests 
that the observational upper limits on ho obtained from 
known isolated pulsars can be used to set a limit on e as- 
suming a value of I zz . However, the moment of inertia is 
very sensitive to the poorly known dense matter equation 



of state (EOS). It can change by a factor of 7 depend- 
ing on the stiffness of the EOS [4(|. Alternatively, with 
Eq. ([2|), one can use Eq. ([I]) to set a limit on the pulsar's 
quadrupole moment without assuming a value of I zz [4~fl ] . 
The limit can in turn set a constraint on the shear modu- 
lus of crystalline color-superconducting quark matter by 
Eq. ((¥]). In particular, with the expression ([5]) for is, we 
can define an exclusion region in the A — fi plane by the 
the following constraint: 
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,(7) 



where ho is the observational upper limit on ho for a 
given pulsar. 

Under the assumption that the pulsar is an isolated 
rigid body and that the observed spin-down of the pulsar 
is due to the loss of rotational kinetic energy as gravita- 
tional radiation, one can also obtain the so-called spin- 
down limit on the gravitational-wave amplitude /i s( j = 
(5GJ zz |/|/2c 3 r 2 /) 1//2 , where /is the time derivative of 
the pulsar's spin frequency [30j| . As it is expected that 
the strain amplitude satisfies ho ^ h s d in general, we 
can derive a constraint on the product A/i based on the 
spin-down limit: 



(A^) sd < 2.1 x 10 5 MeV 2 
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where we have used Eq. (10) of [40j for the moment of 
inertia for strange stars. 

In [3(3] it is reported that the gravitational- wave strain 
upper limit for the Crab pulsar (ho = 3.1 x 10~ 24 ) is the 
closest to the spin-down limit (at a ratio of 2.2). For 
the other pulsars, the direct observational upper limits 
are typically at least 100 times larger than the spin-down 
limits. For pulsars in globular clusters, in which cases the 
spin-down measurement is obscured by the cluster's dy- 
namics, the gravitational-wave observations provide the 
only direct upper limits. In the following, we shall focus 
on the constraints set by the observational data for the 
Crab pulsar. 

The S3/S4 results for the Crab pulsar are plotted in 
Fig- HI In the figure, the solid lines represent Eq. ([7]) when 
the equality holds for the breaking strain er max = 10~ 3 
and 10~ 2 , assuming that M = 1.4M and R = 10 km. 
The pulsar's spin frequency is / = 29.8 Hz and the dis- 
tance is r = 2 kpc. For comparison, the dotted and 
dashed lines are the constraints obtained by the spin- 
down limit [Eq. @], respectively, for cr max = 10~ 3 and 
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FIG. 2: (color online). The dependence of the exclusion region 
on the size of the crystalline core. The mass, radius, and 
breaking strain are fixed at M = IAMq, R = 10 km, and 
o" m ax = 10 -2 . The parameter x = R c /R, where _R C is the 
radius of the core. 



10~ 2 , with the pulsar's spin-down rate |/| = 3.73 x 
1CP 10 Hz s _1 . In the figure, the region to the right of each 
curve is excluded by the corresponding constraint. The 
rectangular box (in dashed line) in Fig.[T]encloses the the- 
oretically allowed ranges of p, and A in compact stars as 
discussed above. It is seen that the direct observational 
upper limits (solid lines) can already place a strong con- 
straint on the possibility that the Crab pulsar might be 
a solid star composed of crystalline color superconduct- 
ing quark matter. For the fiducial values M = 1.4M©, 
R = 10 km, and <7 max = 1CP 3 , the gap parameter A is 
restricted to be less than ~ O(20) MeV for the range of 
quark chemical potential relevant to compact stars. Fur- 
thermore, the more extreme estimate a max = 10~ 2 ex- 
cludes essentially the entire theoretically allowed range 
of A. 

So far we have assumed that the whole star is com- 
posed of crystalline quark matter. However, a more re- 
alistic compact star model would be a crystalline quark- 
matter core surrounded by a nuclear-matter fluid layer, 
which is then followed a thin crust. As a first step to 
investigate the effects of a smaller crystalline core on 
the constraint ([7]), we assume that the quadrupole de- 
formation (and hence the gravitational wave emission) 
is due only to the crystalline core, without consider- 
ing the respond of the outer fluid layer and thin crust. 
The maximum quadrupole moment ((4]) is replaced by 
Qmax ~ l3vcr max R°/GM c , where M c and R c are, re- 
spectively, the mass and radius of the crystalline core. 
To relate Q max to the global parameters M and R of 
the star, it is necessary to know the mass distribution 
of the outer fluid layer, which depends on the nuclear- 
matter EOS. For simplicity, we assume that the density 
falls off as p(r) = p c [l-(r/R) 2 ]/[l-{R c /R) 2 } in the outer- 
fluid layer (i.e., for R c < r < R), where p c is the den- 
sity of the core [42| . The quadrupole moment can then 
be expressed as Q max ~ 13va max R e a(x) /GM , where 
x = R c /R, a(x) = x 6 [l-(l-a; 2 )- 1 (l-3a; 2 /5-2x- 3 /5)], 



and a(l) = 1. Accordingly, the constraint ([7]) is modi- 
fied by multiplying a factor a(x)^ 1 ^ 2 on the right hand 
side of the equation. The results for three different val- 
ues of x are plotted in Fig. H with M = 1.4M , R = 10 
km, and cr max = 10~ 2 being fixed. The figure shows that 
a smaller crystalline core would make the extreme es- 
timate er max = 10~ 2 more compatible with the current 
upper limit for the Crab pulsar. However, the effect due 
to the respond of the outer fluid layer to the internal 
deformation deserves further investigation. 

Summary and discussion. In this paper, we point 
out that the S3/S4 runs of the LIGO and GEO 600 
gravitational- wave detectors are already sensitive enough 
to put an interesting constraint on solid quark-star mod- 
els composed of crystalline color superconducting quark 
matter. We have estimated that the maximum equato- 
rial cllipticity sustainable by these stellar models could 
be as large as ~ 10~ 2 . This is about 2 orders of mag- 
nitude larger than that obtained by Owen for solid 
quark star models composed of quark-clusters, without 
taking into account the effects of color superconductivity. 
We have used the direct observational gravitational-wave 
strain upper limit for the Crab pulsar in the S3/S4 runs to 
constrain the poorly known gap parameter A in the crys- 
talline color-superconducting phase. For cr max = 10~ 3 , 
we conclude that the Crab pulsar could be made entirely 
of crystalline quark matter if A is less than O(20) MeV. 
For the extreme estimate cx max = 10 -2 , the Crab pulsar 
is unlikely to be a complete solid quark star. However, it 
could still contain a smaller crystalline quark core in this 
case. The direct observational upper limit for the Crab 
pulsar should beat the spin-down limit in the fifth science 
run of the LIGO detectors [3Cj • This promises to put a 
stronger constraint on the theoretical models considered 
in this paper. 

Finally, we conclude with a few remarks. (1) It should 
be noted that the existence of crystalline quark matter 
inside compact stars is still a matter of debate. As men- 
tioned before, this phase of quark matter has only been 
studied in phenomenological models of QCD and no anal- 
ysis on the gravitational stability of such stellar configu- 
ration has been carried out. Our work suggests a way to 
put constraint on such exotic compact star models based 
on gravitational-wave observation. We also note that 
Anglani et al. [43| have recently studied the cooling rate 
of these compact stars. This can provide another obser- 
vational constraint. (2) A general argument against the 
existence of quark stars is that (fluid) quark-star mod- 
els are inconsistent to the behavior of pulsar glitches [44[ ■ 
The standard explanation for glitches in the conventional 
neutron star model involves the pinning and unpinning 
of large numbers of superfluid vortices to the solid crust 
[45l l46l |47| . But the significant development of the ef- 
fective theory of QCD over the past few years suggests 
that the quark matter inside compact stars is likely to be 
simultaneously superfluids and rigid solids, which are the 
two crucial conditions for the standard model of glitches 
(see [f| Hll for discussion and [48| for the study of the 
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superfluid mode in the crystalline quark-matter phase). 
(3) One uncertainty in our estimate is the value of the 
breaking strain cr max . One should keep in mind that the 
values <7 max = 10~ 3 — 10~ 2 used in this paper probably lie 
in the upper end of the theoretical range [36], l37f . We also 
assume that the stellar material is everywhere strained to 
the maximum in the analysis. As pointed out by Owen 
[Tj| , since a given pulsar may not be strained to the max- 



imum, thus no upper limit on e max can ever rule out any 
theoretical model. But the credibility of the model will 
face serious challenge as the number of observed pulsars 
increases and tighter limits on e max are placed. 

The author thanks Pui Tang Leung for careful reading 
of the manuscript. This work is supported in part by the 
Hong Kong Research Grants Council (grant no: 401905 
and 401807) and the Chinese University of Hong Kong. 
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